Pressure induced topological quantum phase transition in Sb2Se3 



o 

Oh 
< 

(N 






cd 



c 
o 
i o i 

(N 
> 

m 
m 
m 

o 



x 



Wei Li, 1 - 2 Xin-Yuan Wei, 1 Jian-Xin Zhu, 3 C. S. Ting, 4 and Yan Chen 1 

department of Physics, State Key Laboratory of Surface Physics and 

Laboratory of Advanced Materials, Fudan University, Shanghai 200433, China 

2 State Key Laboratory of Functional Materials for Informatics, 

Shanghai Institute of Microsystem and Information Technology, 

Chinese Academy of Sciences, Shanghai 200050, China 

3 Theoretical Division and Center for Integrated Nanotechnologies, 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

* Texas Center for Superconductivity and Department of Physics, 

University of Houston, Houston, Texas 77204, USA 

(Dated: April 15, 2013) 

Based on the first-principles band structure calculations, we theoretically study the pressure effect 
on a trivial topological insulator (TI) Sb2Se3 and show that there is a topological quantum phase 
transition occurring from a trivial TI to a non-trivial strong TI at a critical pressure value P c ~ 0.35 
GPa. This pressure induced topological quantum phase transition is confirmed by the bulk energy 
gap evolution as a function of the pressure, energy band structure with "band inversion", the Z2 
topological invariant, and the existence of the Dirac-like topological surface states. Our calculations 
can be tested by both spectroscopy and transport experiments. 
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The discovery of time reversal invariant (TRI) topolog- 
ical insulators (TIs) [l|,|2| has triggered intensive interests 
in condensed matter community and materials physics. 
The TIs in two [3, |4| or three (5l4l5J dimensions have 
a bulk energy band gap, like an ordinary insulator, but 
have conducting gapless edge or surface states protected 
by time reversal symmetry. According to the Kramers 
theorem, the edge or surface states, consisting of an odd 
number of massless Dirac cones, are robust against TRI 
perturbations, while the ones consisting of the even num- 
ber of massless Dirac cones are not protected. This even- 
odd effect becomes the key reason why the TI is charac- 
terized by a Zi topological invariant |5j, |l6|. On the one 
hand, when applyi ng the p ressure on the i? 2 non-trivial 
TIs, such as Bi 2 Te 3 17H19J and Bi 2 Se 3 |2y], the supercon- 
ducting phase will emerge. On the other hand, Sb2Se3 is 
known to be a trivial TI in Bi 2 Sc 3 family compounds at 
the ambient pressure. Therefore, there is a question aris- 
ing as to what will happen when applying the pressure 
on a trivial TI Sb 2 Se3. 

Very recently, the pressure induced transition from a 
trivial TI to a non-trivial TI in Sb 2 Se3 was observed by 
the Raman experiment [211 ]. It is shown that this sys- 



number in the two phases divided by the bulk energy gap 
closing point, (iv) The gapless Dirac cone like topological 
surface state inside the bulk energy gap only exists in one 
phase, while the other one does not exist. In addition, 
the mechanism of the pressure induced TQPT observed 
by Raman experiment 
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remains unclear. 

In this Letter, we theoretically study the pressure in- 
duced TQPT in Sb 2 Se3 motivated by aforementioned 
discussions of four fundamental features as well as the 
mechanism for TQPT based on the first-principles calcu- 
lations. Firstly, our results show that when the pressure 
is increased to a critical value P c ~ 0.35 GPa, which is 
much smaller than the value of P^w 2 GPa reported by 



tcm undergoes a topological quantum phase transition 
(TQPT) from a trivial TI phase to a non-trivial one by 
tuning the pressure. From the theoretical point of view, 
to reach such a situation of TQPT, four criteria must be 
satisfied: (i) The bulk energy gap must be closed when 
the TQPT occurs, namely when the pressure is increased 
to a critical value P c , the bulk energy gap should van- 
ish, (ii) The conduction and valence bands cross together 
forming a "band inversion" at the TRI k point in the Bril- 
louin Zone (BZ) at a pressure above the critical value P c . 
(iii) The Z 2 topological invariant should have a different 



Raman experiment 21( , the bulk energy gap is closed ac- 
companied by the TQPT. Secondly, our study also shows 
that a "band inversion" appears at the TRI V point in BZ 
at a pressure value above P c , which is attributed to the 
pressure induced enhancement of the crystal field split- 
ting, resulting in the anti-crossing between conduction 
and valence bands. This is the basic physical picture of 
TQPT induced by the pressure. In order to quantita- 
tively confirm the TQPT indeed occurring at the crit- 
ical pressure value P c ~ 0.35 GPa, we calculate the Z 2 
topological invariant in two phases separated by the bulk 
energy gap closing point and show that the low pres- 
sure phase has a trivial Z 2 topological invariant number 
(0; 000), while the latter one has the non-trivial Z 2 topo- 
logical invariant number (1;000). Finally, according to 
the idea of the bulk-edge correspondence of the TIs, the 
surface spectral function calculations show that the Dirac 
like topological surface states inside the bulk energy gap 
only exists in the high pressure phase, a non-trivial strong 
TI. The test of these calculations should be readily acces- 
sible to both spectroscopy and transport measurements. 
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FIG. 1: (Color online) (a) Crystal structure of the Sb2Se3 
and the primitive translation vectors ti.,2,3 are shown, (b) 
The pressure dependent optimized lattice constants and the 
internal coordinates of the Sb2Se3, where the Se2 site is set to 
be at the origin (0, 0, 0), the two Sb sites are at (±/x, ±/i, ±/x), 
and two Sel sites are at (±i/, ±f, ±f), defined in the unit of 
primitive translation vectors, as shown in (a). 



The Sb2Sc3 has a rhombohcdral crystal structure with 
space group D^ d (R3m), as shown in Fig. Q^a). In Fig. 
[IJa), one can see that the system contains quintuple 
layers ordered in Sel-Sb-Se2-Sb-Sel sequence along the 
c-axis, and the primitive translation vectors ^2,3 read 
as ti = (— a/2, — V3a/6,c/3), t 2 = (a/2, — V3a/6, c/3), 
and t3 = (0, V3a/3, c/3), respectively, where a and c 
are lattice constants in the hexagonal cell. The inter- 
nal coordinates of atoms Sel, Se2, and Sb are set to be 
at the (±1/, ±1/, ±v), (0,0,0), and (±fi,±fi,±fi), respec- 
tively. Without loss of generality, the lattice constants 
and internal coordinates under the pressure ranged from 
GPa to 8 GPa without considering the crystal struc- 
ture phase transition [lj], [l8|, [2l[ occurring are all opti- 
mized based on the first-principles calculations [22j im- 
plemented in the VASP code [24j. The optimized results 
including lattice constants and internal coordinates are 
all shown in Fig. [ljb). At the ambient pressure P= 
GPa the optimized lattice constants for Sb2Se3 are con- 
sistent with those reported by Zhang et al. [9j. From 
Fig. QJb), wc notice that the lattice constants a and c 
are shrunk with the increasing pressure, as one expects 
intuitively. However, it is interesting to note that the 
internal coordinate of Sb atom increases with the higher 
pressure, which is opposite to what we expect intuitively. 
However, it can be understood from the binding energy 
consideration: the coupling is very strong between two 
atomic layers within one quintuple layer (QL), while the 
coupling is much weaker between two QLs predominated 
by the van der Waals type interaction. Consequently, the 
Sb atoms keep the balance position to minimize the total 
energy. 

We now turn to discuss the relationship between the 
bulk energy gap and the pressure. As shown in Fig. [2j 
the calculated bulk energy gap for Sb2Sc3 at the ambient 
pressure is 0.24 eV, which is also in good agreement with 
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FIG. 2: (Color online) The Phase diagram of the bulk energy 
gap as a function of the pressure. The red-dashed line denotes 
the critical point P c ~ 0.35 GPa of the pressure, at which the 
bulk energy gap is closed. The low pressure phase is a trivial 
TI, while the high pressure phase is a nontrivial TI. 



previous experimental and computational results (about 
0.2 - 0.3 eV) H, |H HI. With the increasing pressure, 
the bulk energy gap is gradually reduced. Eventually, the 
bulk energy gap is closed when the pressure is increases 
to a critical value P c = 0.35 GPa. Due to the fact that 
the spin-orbit coupling strength in Sb2Se3 is not strong 
enough to induce the band inversion at the ambient pres- 
sure, Sb2Se3 is maintained as a trivial TI phase at low 
pressure up to P c . If the pressure is further increased, 
the bulk energy gap will be reopened and the system en- 
ters an insulating phase again. The bulk energy gap can 
reach up to a maximal value E g ~ 0.4 eV with the in- 
creasing pressure, which is much larger than the energy 
scale of room temperature. Such a gap closing-reopening 
transition by tuning the pressure indicates that there is 
a phase transition occurring at the critical point P c . 

In order to understand the nature of the new phase as 
aforementioned, wc also calculate the bulk energy band 
structure. The energy band structure with and without 
spin-orbit coupling for the typical pressure values, such 
as at the ambient pressure and the pressure value with 
P= 2 GPa, are all shown in Fig. [3] From Fig. [3]we notice 
that the inclusion of spin-orbit coupling does not lead to 
a band inversion at the ambient pressure. By comparing 
with Fig. [3Jb) and Fig. EJd), one can see clearly that 
the qualitative change is an anti-crossing feature around 
the r point, which indicates an inversion between the 
conduction band and valence band, suggesting that the 
new phase in Sb2Se3 under high pressure is a non-trivial 
TI phase. Therefore, we can conclude that the pressure 
induced phase transition is indeed a TQPT from a trivial 
TI phase to a non-trivial TI phase. 

We now proceed to elaborate the mechanism of the 
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FIG. 3: (Color online) Energy Band structure for Sb2Se3 un- 
der the pressure P= GPa (a)-(b) and P= 2 GPa (c)-(d). 
The spin-orbit coupling is taken into account in (b) and (d), 
while the spin-orbit coupling in (a) and (c) are not included. 
The dashed line indicates the Fermi level. 



pressure induced TQPT. It is known from the atomic 
physics that the spin-orbit coupling strength is propor- 
tional to the atomic number Z 4 (Z = 51 for Sb and 34 
for Se). As we have shown above, the spin-orbit coupling 
strength in Sb2Se3 is not strong enough to induce the 
band inversion at the ambient pressure. Therefore, the 
spin-orbit coupling is unable to induce the band inversion 
even under the high pressure value and it cannot explain 
what the experiment observed [21| either. To explore 
what is the mechanism of TQPT, let us start from the 
crystal filed splitting without considering the spin-orbit 
coupling effect. According to the point-group symmetry, 
the p z orbital is split from the p x and p y orbitals and the 
latter two orbitals remain degenerate on the Sb and Se 
atoms. Due to such splitting, the energy levels around 
the Fermi energy forming the valence and the conduc- 
tion bands exhibit a large energy gap. This gap is much 
larger than the energy scale of the spin-orbit coupling 
at the ambient pressure (compare Fig. [3ta) with Fig. 
EJb)). When the pressure is applied, the crystal structure 
will be squeezed. As a result, the level splitting will be 
enhanced and the bulk energy gap between the valence 
band top and the conduction band bottom at T point 
will be gradually reduced. Eventually, the valence band 
and the conduction band may anti-cross to form a semi- 
metallic state when the pressure is increased to a large 
value, see Fig. |3jc). In Fig. Etc), we notice that there are 
two accidental degenerate states appearing around the T 
point, which is not protected by time reversal symmetry. 
Consequently, if we further take the spin-orbit interac- 
tion into account, these accidental degenerate states will 
be broken and a non-trivial TI is formed (see Fig. GD^d)). 
To firmly confirm the TQPT nature induced by prcs- 



TABLE I: Parities 8 t at eight TRI k points for Sb 2 Se 3 under 
different pressure values. The Z2 topological invariants are 
(1;000) for a pressure above the critical value P c while the 
Zi topological invariants are (0; 000) for a pressure below the 
critical value P c , which indicate a non-trivial strong TI and a 
trivial TI, respectively. 

8\ 8 2 S 3 84 8 5 5e 87 8$ (vo\viV2V'a) 



OGPa +1 +1 +1 +1 

0.2GPa +1 +1 +1 +1 

0.3GPa +1 +1 +1 +1 

0.4GPa -1 +1 +1 +1 

2GPa -1 +1 +1 +1 

5GPa -1 +1 +1 +1 



+ 1 +1 +1 +1 (0;000) 

+ 1 +1 +1 +1 (0;000) 

+ 1 +1 +1 +1 (0;000) 

+ 1 +1 +1 +1 (1;000) 

+ 1 +1 +1 +1 (1;000) 

+ 1 +1 +1 +1 (1;000) 



sure, we further calculate the Zi topological invariant. 
Because the existence of the spatial inversion symmetry 
in Sb2Se3, the calculation of the Z^ topological invariant 
can be dramatically simplified by the so-called "parity 
method" [2jj. Accordingly, the Zi topological invariant 
of Sb2Se3 under the pressure environment can be ob- 
tained from the wave function parities at the eight TRI k 
points defined as r i= ( ni „ 2 „ 3 ) = \{n\Qt\ + niGi+nzGz), 
where Gj are primitive reciprocal lattice vectors with 
11 j = or 1, and determined by the quantities 



N„ 



I &m (Ti). 



(1) 



Here, £ 2m (r,) = (<5>im,r z \P\<5>2 m ,r.) = ±1 is the eigen- 
value of parity operator P corresponding to even (odd) 
parity of the Bloch functions at the 2m-th occupied band 
and the TRI k point Ti (i = 1, ..,8). The sum is taken 
over all of the occupied bands N occ with only even band 
index due to the Kramers degeneracy ^ m = ^2 m _i at 
the TRI k point. In Sb2Se3 system, there are four inde- 
pendent invariants (vq; V1V1V3), 
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n k = l,rij^ k =0,l 



where vq is the strong topological invariant and indepen- 
dent of the choice of primitive reciprocal lattice vectors 
Gj while V\ , Vi , and ^3 defined as weak topological invari- 
ants are not. The combination of these four independent 
invariants (i/q; v\vivz) clearly distinguish three classes of 
states: non-trivial TI with strong TI and weak TI and 
trivial TI. The calculated parities Si (i = 1, ...,8) at the 
eight TRI k points with relative coordinates (0,0,0), 
(0,0,0.5), (0,0.5,0), (0,0.5,0.5), (0.5,0,0), (0.5,0,0.5), 
(0.5,0.5,0), (0.5,0.5,0.5), respectively, in a primitive 
reciprocal-lattice are listed in Table HI The Zi topologi- 
cal invariant (i/q\ vivivz) can be obtained by multiplying 
over the parities Si of all TRI k points according to Eq. 




FIG. 4: (Color online) Energy and momentum dependence of 
the surface spectral function for Sb2Se3 with ambient pressure 
(a) and pressure value P= 2 GPa (b) on the [111] surface. The 
red regions indicate bulk energy gaps while the green regions 
indicate bulk energy bands. The Dirac-like topological surface 
states can only be seen clearly around the V point under the 
pressure P= 2 GPa (lower panel). The dashed line indicates 
the Fermi level. 



([2]) and are also listed in Table Q] One can see that the Z-i 
topological invariants are (1;000) for Sb2Se3 at a pres- 
sure above the critical value P c indicating a non-trivial 
strong TI while Z2 topological invariants arc (0; 000) for 
Sb2Sc3 at a pressure below the critical value P c indicat- 
ing a trivial TI. This agrees with the discussions above 
that a "band inversion" between the conduction band 
and valence band occurs at the T point. It is also shown 
that there is a TQPT occurring at the critical pressure 
value P c = 0.35 GPa. 

Finally we discuss the topological surface states based 
on the idea of the bulk-edge correspondence of the 
TIs [28[. In order to reveal the topological features of 
the Sb2Sc3 system under the pressure, wc study the topo- 
logical surface states of a semi-infinite system with and 
without the pressure based on the surface Green func- 
tion [29j,|30[ constructed by the maximally localized Wan- 
nier function fitting to the first principles band structure 
calculations using the method developed by Marzari and 
co-workers [3l|, |32|. The dispersion of the surface states, 
as shown in Fig. 21 are obtained by evaluating the imag- 
inary part of the surface Green function. By comparison 
of Fig. HJa) with Fig. IDJb), one can clearly see the topo- 
logical surface states, which form a single Dirac cone at 
r point, only exist at the pressure P= 2.0 GPa. This in- 
dicates a non-trivial TI nature and agrees well with the 
bulk parity analysis and band structures calculations. 

In conclusion, we have performed the first-principles 
calculations for the pressure effects on Sb2Sc3 and shown 



that there is a TQPT occurring at a critical value P c ~ 
0.35 GPa. To firmly confirm such TQPT, wc firstly cal- 
culate the bulk energy gap evolution as a function of the 
pressure and show the bulk energy gap closc-rcopening 
transition. Secondly, the band structure calculations 
show that there is a "band inversion" appearing at the T 
point in BZ at a pressure above P c , which is attributed to 
the pressure induced enhancement of crystal field split- 
ting, resulting in the anti-crossing of the conduction and 
valence bands. Thirdly, the calculated Z-i topological in- 
variant indicates that a trivial TI phase (0; 000) and non- 
trivial strong TI phase (1; 000) is separated by the critical 
pressure P c . Finally, the calculated Dirac-like topologi- 
cal surface states only exist in the non-trivial strong TI 
phase (1; 000) region. These results can be tested directly 
by spectroscopy techniques, such as ARPES and STM, 
and transport measurements. 
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